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We report on numerical simulations of one dimensional maximally supersymmetric SU(N) Yang- 
Mills theory, by using the lattice action with two exact supercharges. Based on the gauge/gravity 
duality, the gauge theory corresponds to N DO-branes system in type IIA superstring theory at 
finite temperature. We aim to verify the gauge/gravity duality numerically by comparing our 
results of the gauge side with analytic solutions of the gravity side. First of all, by examining the 
supersymmetric Ward-Takahashi relation, we show that supersymmetry breaking effects from the 
cut-off vanish in the continuum limit and our lattice theory has the desired continuum limit. Then, 
we find that, at low temperature, the black hole internal energy obtained from our data is close to 
the analytic solution of the gravity side. It suggests the validity of the duality. 
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1. Introduction 

Gauge/gravity duality asserts an equivalence between strongly coupled gauge theory and the 
classical gravity on curved space, which was originally stated as AdS/CFT correspondence which 
includes the supersymmetry by Maldacena[jl]]. From the duality, we expect that strongly coupled 
gauge theories, which are usually difficult to calculate by hand, can be analytically solved via the 
gravity side. So, there are many applications from the context, getting over the barrier among 
fields (for example, elementary particle physics, cosmology, condensed matter physics and so on). 
However, it is a conjecture and therefore verifying the duality in some way is desirable. 

We aim to verify the gauge/gravity duality from lattice simulations in one dimensional super- 
symmetric Yang-Mills theory with sixteen supercharges. The theory is obtained by dimensional 
reduction from lOd jV = 1 SYM (or 4d JV = 4 SYM). Actually, we introduce temperature into 
the theory by imposing the anti-periodic boundary conditions on fermions. Based on the duality, 
the gauge theory corresponds to N DO-branes system in type IIA superstring theory at finite tem- 
perature. In particular, at low temperature, the gauge theory becomes strong coupling and using 
analytic techniques to examine the dual black hole physics from gauge theory become difficult. 
So, we use the lattice gauge theory to analize the gauge theory. From comparisons between lat- 
tice results and analytic solutions of the gravity side, we discuss the validity of the gauge/gravity 
duality. 

There are two previous works about numerical simulations of the Id maximally supersym- 
metric Yang-Mills theory: non-lattice simulations done by Nishimura et al.^, ||, |], ||] and lattice 
simulations done by Catterall and Wiseman^]. Supersymmetry is broken by the cut-off effects 
in their regularized theories. Nevertheless, both results are consistent with the gravity side from 
UV-finiteness of Id gauge theories. In contrast, we employ the lattice formulations with a few 
supersymmetric charges on the lattice [[7], |h, which have been recently developing, in particular, 
our lattice theory has two exact supercharges even on the lattice. We expect that the lattice theory 
has some advantages, for example, clear signals etc., thanks to the exact charges, in high accuracy 
verifications of the duality. 

In section 2 we explain our lattice theory and then in section 3 we see some details of simu- 
lation techniques. In section 4 we show that the lattice theory has the correct continuum limit by 
computing the supersymmetric Ward-Takahashi relation. In section 5 we show the internal energy 
of the dual black hole obtained from gauge theory side as an evidence of the gauge/gravity duality. 

2. Id SYM with 16 supercharges 

The supersymmetric SU (N) Yang-Mills theory with sixteen supercharges is a gauge theory in 
which a gauge field of the temporal direction A\ interacts with nine scalar fields Xj(i = 1, • • • ,9) 
and real sixteen fermions yf a (a = !,•■■ , 16). The continuum action is given by 



cont 




(2.1) 
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where A is the 't Hooft coupling constant. Here, all fields are expanded as (p(t) = L^Lj 1 (p a {t)T a 
by gauge group generators T a of the SU (N) group 1 . Also, the covariant differential operator D\ 
is defined through D\ <p = d\(p + i\A\ , <p]. The y, are real symmetric matrices which satisfy the nine 
dimensional Euclidean Clifford algebra. 

The realization of supersymmetry on the lattice has been a difficult issue due to the lack of 
Leibniz rule on the lattice for a long time. However, recently, Sugino proposed a lattice formu- 
lation of maximally supersymmetric Yang-Mills theory with two exact supercharges[[7]] from the 
topological twisted version In the twisted theory, the original action eq.(ll) can be rewritten 
as a closed form, S cont = <2+<2-(---)> using two supercharges Q± where <2± are gauge transfor- 
mations. From the nilpotency of Q± up to gauge transformations, we can see that the action is 
invariant under (^-transformations, without the obvious use of Leibniz rule. 

Let us consider one dimensional lattice of the size L with the periodic boundary condition. 
Scalars and fermions are defined on sites labeled by t = 0, •■ • ,L — 1 while a gauge field is defined 
on links through the link field U\ to realize the exact gauge invariance. Our lattice action is defined 
by 



N 



L-l 



-if>(0Si(0-5 I ei]kBi(t)\Bj(t),B k (t)] 
i=\ J ;j,*=i 



4 3 



!=1 



(2.2) 



where Ao is a dimensionless 't Hooft coupling constant defined by Ao = Xo? with the lattice spacing 
a. Here, B,-,C, 0± and y±n,y]±,E,±i are some combinations of original scalars and those of fermions, 
respectively (our notation follows [Q], or see [13 ]). Lattice Q± -transformations are defined over the 
new variables 2 , 



Q±U 1 (t) = ixi/ ±1 (t)U 1 (t), Q±V±i(t) = iy±i(t)y±i(t)±iV 1 <l>±(t), ••• . 
From the definitions above, Q 2 ± are gauge transformations with gauge parameters (j)± , C, 



Q\ = , Q 2 = 8m , {Q + ,Q-} = 5_,c, 



(2.3) 



(2.4) 



where 8 m is a gauge transformation with the parameter CO. As a result, Q± -invariance is realized 
even on the lattice, because of eq.(2.2) and the exact gauge invariance of the lattice theory. 

The continuum limit is realized by taking Xq to zero while keeping a typical scale of this 
system (e.g., the dimensionful 't Hooft coupling A). Also, the lattice action has no doublers because 
(DqD q )(p) = 4sin 2 (/?i/2) where D Q is the free limit of the lattice Dirac operator. 

To introduce temperature, we change the boundary condition on fermions from periodic one 
to anti-periodic one, while keeping that on bosons, 

U l (t) = U l (t + L), X i (t)=X i (t+L), (i=l,-,9), (2.5) 
Wa(t) = -Ya(t + L), (a = l,---,16), (2.6) 

where temperature is defined by T = 1 /(La). Temperature breaks all supersymmetries explicitly 
and, of course, <2±-invariance. Hereafter, T denotes a dimensionless temperature T/A 1//3 . 

'The generators T a satisfy the normalization condition, tr(T a T b ) = 8 a t,. 
2 The transformations here are those for U\ and \ff±i (for the others, see [|13|] ). 
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3. Simulation details 

We use the standard Hybrid Monte Carlo method. But, there are two additional difficulties in 
our fermion sector: the 4-fermi interaction and the pfaffian, as explained below. 

After Q± -transformations, the action eq.(Oh includes a cut-off order 4-fermi interaction as 



(2 + <2_tr(y/+iV-i) ~tr ({!//•+!, i//--i} 2 ) — ► tr (o 2 + o{\i/+ 1, Y-i}) ■ (3-1) 



For the 4-fermi interaction, we introduce an auxiliary field a to write it as the third term in eq.( 3. 1 ). 
We treat tr(a 2 ) as a part of the boson action and tr(a{i// +1 , V-i}) as a part of the fermion bilinear 
Sp = Y,t Y T W^VtO' without integrating a. So, we regard 10+1 bosonic fields, U\,A%, • ■ ■ ,Aio, d, 
as configurations generated by HMC method. 

The integration of fermions becomes the pfaffian, pf(£>), which generally takes complex val- 
ues. We treat the absolute value and the complex phase of the pfaffian, individually, to avoid 
the sign problem. The absolute value can be given as an integral by pseudo fermion since 
|pf(D)| = det(D t D) 1 / 4 , and the 4th root is approximately by the rational expansion, 



|pf(D)|= /D^expl-^VW 



■ 



where the order M and the coefficients of the approximation are determined from the range 
of D f D's eigenvalues measured in the simulations. We compute inversions of D'D with shifts /3, 



in eq.(p.2[) by using multi-mass shifted solver. Also, for the phase of pfaffian, we use the phase 
quench, or use the phase reweighting method if we want to include the effect in the results. 

For r> 1, HMC method stably works and we can obtain sufficient statistics. However, as 
temperature decreases, at some temperature (which actually depend on AO, the magnitudes of scalar 
fields monotonically increase against Monte Carlo trajectories and therefore the fhermalization 
does not occur. This instability is related to the classical flat direction of the boson action. To 
avoid the instability, in parameter regions where run-away modes to the flat direction appear, we 
introduce a mass term of scalar fields, 

Smass = Mo £ E tr , (3-3) 

where jUo is a dimensionless mass. Hereafter /i denotes a dimensionless mass /i/A 1//3 . 

4. Supersymmetric Ward-Takahashi relation 

Supersymmetry is broken by the lattice cut-off. In the classical continuum limit, the break- 
ing effects by the cut-off identically vanish. By contrast, in the quantum theory, it is not so clear 
whether the breaking effects generally vanish in the continuum limit because of ultra violet diver- 
gences and non-perturbative effects. Fortunately, the Id gauge theory is UV-finite and any opera- 
tors which break supersymmetry are not generated radiatively, so our lattice theory has the correct 
continuum limit, at least, in the perturbation theory. However, we do not know a-priori whether 
the similar argument is possible beyond the perturbation theory. Also, we must divide the cut-off 



effects from the other breaking sources, the temperature and the mass term, eq.(3.3). 
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Figure 1: The ratio < <y[i)y$> * s Pl ottec l against the temporal direction for N = 3 (Left) and N = 4 (Right) 
on the L = 16 lattice. The labels si, • • •, sl6 denote the spin indices of supercurrent J\. Red lines represent 
the fitted masses obtained by fitting plateaus. They are consistent with the mass, pi 2 (black lines) within 
statistical errors. Deviations from the plateaus near the lattice boundary correspond to the contact terms. 



For the issue, Kanamori and Suzuki used a method which can extract only the cut-off effects 
in lattice simulations of 2d JV = (2,2) SU(2) SYM[10]. The method is a simple one measuring of 
the partially conserved SUSY current on the lattice. We use the same method and check whether 
the cut-off effects vanish in the continuum limit without relying on the perturbation theory. 

For the continuum theory with mass term, we have the partially conserved supersymmetric 
Ward-Takahashi relation, 



{diJ la (t)0(s)) = li 2 {Y a {t)0{ S ))-8{t-s)(Q a 0{ S )), 



(4.1) 



where G is an arbitrary operator and Q a {a = 1 , — , 16) are supercharges which generate super- 
symmetric transformations. Here J\ and Y are the supercurrent and the breaking term from mass 
term, respectively, which are defined as 



J - N 



16 



1 



16 



i=l 0=1 
N 9 16 



ij=l 0=1 



(4.2) 
(4.3) 



i=10=l 



In the continuum theory, the supersymmetric Ward-Takahashi relation, eq.(fO|), holds even at 
finite temperature. It means that we can find the cut-off effects by measuring a lattice counterpart 
of the relation. Actually, we compute the following ratio, 



< VjJU^YpiO) > 
<Y a (t)Y p (0)> ' 



for a,j8 = !,-•■, 16, 



(4.4) 



where Vj is the symmetric covariant difference operator. Here we only use the forward covariant 
difference operator in the lattice definition of the supercurrent J\. Also, for fixed a, (5 is uniquely 
determined because correlators in the denominator with other j3 are nearly zero, that is, the ratio is 
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Figure 2: The internal energy of black hole, normalized by A^ 2 . Data points(red circles) are our results for 
N =\A. The dashed green line corresponds to the result obtained by high temperature expansion at next to 
leading order in the large N limit. The dashed blue curve represents the analytic solution of the gravity side. 



meaningless. If the cut-off effects vanish in the continuum limit, the ratio must be \i 2 in the limit, 
from eq.(fll|). 

In Figure [II we plot the ratio for N = 3 and N = 4 with T = 1, /x 2 = 0.01 and the lattice size 
L = 16. The horizontal axis represents the temporal direction. The corresponding lattice spacing 
is a = 0.0625 in the unit of A = 1. For both cases, clear plateaus are observed within statistical 
errors. We perform the constant fit for the plateaus. The obtained values are consistent with \i 2 = 
0.01 within the statistical errors. The result suggests that cut-off effects mostly vanish near the 
continuum limit and our lattice theory has the correct continuum limit beyond the perturbation 
theory for, at least, T = 1 of N = 3 and N = 4. 



5. Internal energy 

The internal energy of the black hole, associated with the black hole thermodynamics, is one 
of simple examples to test the gauge/gravity duality. In the gravity side, the internal energy E is 



related to temperature T through an analytic formula[|l 1|], 



1 { E \ ( T \ 14/5 9 l\ 13 ^/rc\ i4,,1/5 



15 z - } =7.41--- . (5.1) 



N 2 \V/ 3 J U 1/3 7 14 1 W 

We compute E = — ^-lnZ, where j3 is the inverse temperature, from our data and compare it with 
the above analytic formula. 

In Figure ||, we show the internal energy versus temperature for N = 14 and \i 2 = 0. We used 
two different lattice sizes, L = 8 for T > 1 and L = 16 for T < 1. For high temperature, we see 



that our data and the result obtained by high temperature expansion at next to leading order[12] 
agree as expected. However, as temperature decreases, our data departs from the curve of the high 
temperature expansion around T = 1.5 and it is smoothly close to the analytic curve of the gravity 
side, eq.(pTT|). This suggests the validity of the gauge/gravity duality in this system. 
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In N = 14, HMC-runs are stable for T > 0.5 without the mass term, however, for T < 0.5, the 
instability from the flat direction occurs. So, we must use N > 14 to explore possible further lower 
temperature. Also, simulations with different lattice spacings at same temperature are necessary to 
take the continuum limit. Such simulations to obtain high accuracy results are in progress[14]. 
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